Subdiffusive fluctuations of "pulled" fronts with multiplicative noise 
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We study the propagation of a "pulled" front with multiplicative noise that is created by a local 
perturbation of an unstable state. Unlike a front propagating into a metastable state, where a 
separation of time scales for sufficiently large t creates a diffusive wandering of the front position 
about its mean, we predict that for so-called pulled fronts, the fluctuations are subdiffusive with root 
mean square wandering A(t) ~ t 1 ^ 4 , not t 1 ^ 2 . The subdiffusive behavior is confirmed by numerical 
simulations: For t < 600, these yield an effective exponent slightly larger than 1/4. 
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Since the late 1930s, when the concept of front prop- 
agation emerged in the field of population dynamics 
jl],|), interest in this type of problems has been growing 
steadily in chemistry || , physics [0 and mathematics || . 
In physics, the importance of the problem has become 
more and more clear since it plays a role in a large vari- 
ety of situations, ranging from reaction-diffusion systems 
to pattern forming systems in general 0. 

Front propagation into unstable states is an interesting 
dynamical problem by itself. For a front evolving from 
a local perturbation there are but two possible propa- 
gation mechanisms that are determined by the nonlin- 
earities in the equation of motion: Either the nonlin- 
earities determine the velocity of the front that then is 
called "pushed" ; or the nonlinearities simply cause satu- 
ration and the velocity is determined by a linearisation 
about the unstable state. Fronts of this type are called 
"pulled" because they are "pulled along" by the spread- 
ing and growth of small perturbations about the unsta- 
ble state j?]]. Hence pulled front propagation can occur 
only if the penetrated state is linearly unstable. The 
pushed and pulled regimes are also known as nonlinear 
and linear marginal stability S. For the discussion be- 
low, it is important to realize that pushed fronts relax 
exponentially in time to their long time asymptotes, but 
that pulled fronts relax algebraically without characteris- 
tic time scale . Hence an adiabatic decoupling of some 
outer dynamics from the internal relaxation of a pulled 
front is not possible |9fl , and stochastic pulled fronts may 



show anomalous scaling |10|. 

Generally, noise can affect the phenomenological de- 
scription of a reaction-diffusion system in various ways. 
A first possibility is intrinsic noise modelled typically by 
additive thermal noise in a Langevin type equation. A 
second possibility, on which the present paper is focused, 
is at the external level, e.g. due to fluctuations of some 
control parameter. An example are the fluctuations of 
the luminosity intensity in the photosensitive Belousov- 
Zhabotinsky reaction [jll). Such fluctuations enter the 
dynamical equation as multiplicative noise. 

The multiplicative noise of the control parameter usu- 



ally results in a modification of the mean propagation 
velocity of the front and in a stochastic wandering of the 
front position around its mean propagation. This means 
that the noisy front can be thought of as a coherent struc- 
ture whose motion can be decomposed into drift plus 
Brownian motion, very much like a particle sedimcnting 
in a fluid. The drift component corresponds to an aver- 
age front, with the average taken over the ensemble of 
all the realizations of the noise. It propagates according 
to a deterministic equation of motion, whose dynamical 
parameters are in the simplest case just renormalised by 
the noise. Theoretically, the important question then 
arises whether the effects of the fluctuations of the front 
can be understood in terms of a diffusive or subdiffusive 
wandering of some suitably defined front position. 

The renormalisation of the front velocity has been 
studied in the pushed and pulled regime Il3], while the 
wandering process is understood only in the pushed case 
|]i~3| , where it has been shown to be diffusive: the root 
mean square position of the front A grows with time as 
yf2Dft. Actually, the expression for the effective front 
diffusion coefficient D / derived by Armero et al. was 
found to break down for pulled fronts, and it was sug- 
gested that the wandering of pulled fronts is subdiffusive. 

In this paper we take up the issue of the stochastic 
wandering of pulled fronts about their mean position, 
and predict that in the presence of multiplicative noise 
pulled fronts behave subdiffusively, with A ~ i 1 / 4 . This 
prediction is based on two different arguments. First of 
all, we heuristically insert the leading edge asymptotics 
of the relaxing pulled front into the expression for the dif- 
fusion coefficient Df of pushed fronts, and immediately 
find A ~ i 1 / 4 . Our second argument for the subdiffusive 
A ~ i 1 / 4 behaviour comes from mapping the dynami- 
cally important region onto the KPZ equation. We finally 
also present data of extensive numerical simulations that 
support our analytical prediction that the wandering is 
subdiffusive with exponent close to 1/4. 

The qualitative difference between pushed and pulled 
fronts results from the fact that the dynamically impor- 
tant region for pushed fronts is the interior front region, 
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whose extent is finite, while that of pulled fronts is the 
leading edge ahead of the front Q . Starting from a local 
initial perturbation, the leading edge region grows with- 
out bound and as we shall see, this causes the subdiffu- 
sive behaviour. The power law relaxation of deterministic 
pulled fronts is another manifestation of the leading edge 
dominated dynamics of pulled fronts . 

For concreteness, we derive our results by including 
noise in the one dimensional prototype front equation 



/(</>)= 0(1 -0)(a + 0). (1) 



Here a is a parameter which plays the role of the control 
parameter. Equation ([!]) has a stable state 0=1 and 
a stationary state = whose relative stability can be 
tuned by changing the value of the parameter a. The case 
— ^ < a < \ leads to pushed dynamics, while i < a < 1 
produces pulled fronts 0. For the case a — 1, which we 
will study, the so-called Fisher-Kolmogoroff-Petrovsky- 
Piscounoff (F-KPP) equation jl]|2| is recovered. 

Let us assume now that the parameter a is replaced 
by a new fluctuating parameter a(x, t) with average a, 
a — > a(x, i) = a + fi(x,t), where n(x,i) is a Gaussian 
noise with the moments: 

(MOM)>„ = 0, (2) 
(MOr, t) fi(x\ <'))„ = 2e C(X„ \x - x'\) S(t - t'), (3) 

with J dx C(Ap, \x\) = 1. We interpret the stochastic 
p.d.e. defined by (|l|) - (||) in the Stratonovich sense [fTi) . 
Notice that if 1/ A M is much smaller than any other length 
scale in the system, the noise defined by the correlator 
(|J) is effectively white in both time and space. 

Since according to (|l|) cj> converges to 1 and is noise- 
less behind the front, we can suitably define the position 
Xf{t) of a noisy front propagating to the right into the 
unstable state (j> = by 



*/(*) 



dx <p(x, t), 



(4) 



The displacement Ax/(i) = Xf(t) — x/(0) on average 
grows with the noise renormalized mean velocity vr = 
(x/)^. The fluctuations about the mean displacement 
(Ax/(t)) AI = vrX are measured by 



A(t) = y((A X/ (t)-<A X/ (<)),) 2 )^ 



(•5) 



If we relate A(t) to a diffusion coefficient Df by writing 



A 2 (t) = / dt' 2D f (t'). 



(G) 



then for pushed fronts the following expression for the 
diffusion coefficient Df can be derived 0,|l5f: 



D f = e 



(7) 



In this formula, cf> is the deterministic field associated 
with the front moving with the renormalised pushed ve- 



is the derivative of the reac- 

and 



locity v R , g(<p) = 

tion term with respect to the control parameter 
£ = x — VRt is the comoving coordinate. 

For pushed fronts, Df given by (0) is finite and time- 
independent, and hence this gives the diffusive behavior 
A 2 (£) = 2D ft. This means that on sufficiently long time 
scales the random displacement is approximately Marko- 
vian, i.e., the sum of uncorrelated and equally distributed 
random displacements on shorter time scales. 

As an example of a pulled front with multiplicative 
noise, we now study the case a = 1: 



dt dx 2 



(8) 



The noise renormalized mean velocity v* R of the pulled 
front can be calculated explicitly n% : 



v*r = (*(*)>„ = 2^(1+ eC(0)). 



(9) 



However, it is immediately clear that the fluctuation for- 
mula (0) cannot naively be extended to the pulled regime. 

First of all, for a pulled front the expression (^) simply 
diverges. The divergence of solvability-type expressions 
actually holds more generally for perturbative expansions 
about a pulled front ||. For a pulled front, the dynami- 
cally important region is the leading edge defined as the 
region where linearisation about the unstable state is a 
valid approximation; the fact that solvability-type inte- 
grals like (0) diverge there reflects that the dynamically 
important region becomes semi-infinite. 

Second, a pulled front has no characteristic relaxation 
time , so there is no reason for the Markovian approxi- 
mation underlying diffusive wandering. Rather the lead- 
ing edge relaxes asymptotically as @] 

w a£ R e' x *^ R e-^'^/t 5 /' 2 , \* R = v R /2, (10) 



for £r 



v* R t > 1 and t > 1. 



The presence of the a^R/t 3 ^ 2 term in front of the ex- 
ponentials is actually the fingerprint of the full equation 
being nonlinear. The expression (JlCj) defines a time de- 
pendent Gaussian cutoff £ c ~ \/ADt, which regularizes 
the integrals in (Q). In fact, the evaluation of (Q) with 
© yields 



Df(t) 




(t » 1). 



(11) 



Notice that for large times Df(t) vanishes, marking the 
nondiffusive wandering of pulled fronts. Insertion into 
(|) yields 
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A(i) 



df D f {f) 



12s 



1/2 



,1/4 



(12) 



so the fluctuations are subdiffusive with exponent 1/4 
rather than 1/2. 

Although the above argument does capture the essen- 
tial features of fluctuating pulled fronts, it is not entirely 
systematic, as it is based on the extrapolation of the solv- 
ability condition (0) to the pulled regime. 

In order to substantiate the scaling A(t) ~ i 1 / 4 for 
a relaxing pulled front with a time-dependent analysis, 
let's go back to Eq. (||). The leading edge region can be 
studied by means of the leading edge transformation, 



£ = x - v*t , v* = 2 
Eq. (§h can then be written as 

at =D W~ i ' 



(13) 



A* = 1. 



(14) 



+ e* [(1 + fi)i'e-t - fitl) 2 e- 2i - V 3 e~ 341 
For £ ^> 1, the nonlincaritics can be neglected 



dip _ n d 2 tp 

~dt = 



D^j+fiip, for 



(15) 



Notice that the noise in this "directed polymer" equation 
still is multiplicative. The Cole-Hopf transformation 



converts (|l5|) into an equation with additive noise: 



dh 
~dt 



D 



d 2 h 



D 



dh 



H, for £>1. 



(16) 



(17) 



Eq. ( |17| ) is the celebrated 1-dimcnsional Kardar Parisi 
Zhang (KPZ) interface equation |16)| . 

The essential difference between our problem and pre- 
vious studies of the KPZ equation are the initial and 
boundary conditions. After some temporal evolution, the 
nonlinearities in the original <f> equation will lead to the 
fluctuationless saturation of <p a t the value of unity for 
£ <C — 1, which corresponds to the fluctationless slope 
h w A*£ behind the front: It is as if the KPZ equation 
has to be solved in the positive half-space with (roughly) 
a fixed boundary. On the other hand, by translating jTo| ) 
back into h, we see that for large £ and t, the average in- 
terface shape h av should be given by 



ln(«W* 3/2 ) + A*£ - \* R £ R - &/4Dt 



(18) 



Thus, apart from the logarithmic term the average inter- 
face is essentially tilted but flat up to the time-dependent 
cross-over £ c s» y/4Dt |l^] , and beyond £ c it has the shape 
of a downward curved parabola with time dependent cur- 
vature. Together with the fact that the nonlinear term 



in (|17j) gives an average nonzero growth velocity, this 
makes the problem into a nonstandard fluctuating in- 
terface problem. Our central approximation is now to 
consider the relaxing front in the essentially straight but 
fluctuating section between and v 4Dt as a KPZ in- 
terface with time dependent length L = 0(£ c ). As the 
scaling exponents of the KPZ equation are robust with 
respect to a geometric change of the fluctuating surface 
p8[ , we use the KPZ scaling functions for the root mean 
square width W of the interface h, 



W(L,t) = t^Y(±y 13=1/3, 



z = 3/2, (19) 



with the bar denot- 



where W = \J (h(x,t) - h{x,t)) 
ing a spatial average. The scaling function Y(s) will de- 
pend on the shape of the roughening surface, but always 
has the limits Y(s) — > s^ 13 for s — * oo, Y(0) as const. 
Inserting our approximation L ~ we get: 



W(L,t) ~L zf} - (Vt) z/3 = i 1/4 . 



(20) 



The final step of our argument is to convert this result in 
a prediction for the fluctuations of the front position. If 
we measure the position of the front by tracking a certain 
height c, cj)(x c ,t) = const = c, and use the relations jl^ ) 
and (U), we find: 



(21) 



This implies that fluctuations in h are just identical with 
fluctuations in x c . Therefore we get 



(22) 



A(t) - t 1/4 



which reproduces the scaling of our previous result (112 

We have also performed numerical simulations of the 
noisy front equation (|l|) with a = —0.3 (pushed) and 
a = 1 (pulled, F-KPP Equation (||)) following the lines 



of 1 13 1 . The initial condition was taken as a step function 
<p(x,0) = 0(xq — x). The numerical integration has been 
performed using a standard explicit Euler algorithm, in 
both cases the value of the noise was set to e = 0.5, and 
the zero value of the spatial noise correlator C(0) was cho- 
sen as the inverse spatial integration mesh, C(0) = 1/ Ax 
Jl3| . The result is shown in Fig. 1, where the function 
A(f) is plotted in both the pushed and the pulled case. 

The specific features of the pulled regime make the 
problem quite delicate from the numerical point of view. 
In order to minimize finite size effects, which are partic- 
ularly worrisome in this regime Q , we have worked with 
a large system size (L — 3000) and gridsize Ax — 1 (the 
change in v* and D due to the finite gridsize effect was 
taken into account following the prescription of fjj ) . This 
made sure that even at time t — 600, the leading edge of 
the front never reached the boundary of the system. 

We have also checked our program and system size ex- 
tensively both for deterministic and noisy fronts, taking 
into account grid and time step effects according to M. 
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ponent due to these special initial conditions is an issue 
that will be addressed elsewhere. 

We finally stress, that our results apply to a much 
larger class of equations than nonlinear diffusion equa- 
tions (|l|). The methods of generalization are analogous 
to those of a closely related result is the general 

argument put forward in fTcfl that noisy pulled fronts in 
more than one dimension should not obey KPZ scaling. 
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FIG. 1. Diffusive and subdiffusive spreading of the front po- 
sition. The dot-dashed curve correponds to the pushed case 
(o = —0.3) and the solid one corresponds to the pulled case 
(a — 1). The dashed straight line is the prediction dl2]), while 
the dotted line indicates a slope 1/2. 

Our final result, based on averaging over 10000 front 
realizations, is shown in Fig. 1; it clearly confirms the 
subdiffusive behaviour predicted by our analytical argu- 
ments. Quantitatively, when we associate a single ef- 
fective exponent with the late time slope in the log-log 
plot of Fig. 1, we get an effective exponent of about 0.29 
rather then 1/4. Over the time interval we have stud- 
ied, the actual value of A(t) is somewhat larger than 
an asymptotic prediction (|12|), which is indicated with 
a dashed line. This may be due to the fact that ( |l2] ) 
only gives the behavior for such long times that the time 
integral is dominated by its large t behavior. The fact 
that A is only of the order of 4 at our latest times sug- 
gests that this asymptotic regime is only reached at very 
late times. Indeed, assuming that finite size effects are 
negligible, we attribute the fact that the effective expo- 
nent is slightly larger than 1/4 to the presence of slow 
crossovers, which surely are present in the system. Some 
of these can be estimated, while others are more difficult 
to trace, (i) We already noticed previously that we are 
actually dealing with a slightly curved KPZ interface, for 
which the crossover scaling functions are not known, and 
that the way in which the cutoff £ c = 0{^/t) enters the 
KPZ analysis requires further study, (ii) The corrections 
to our asymptotic estimates for the integrals in (0) are all 
of order 1 / \ft, with possible logarithmic corrections 
This indicates that the corrections to the scaling A ~ i 1 
are of order i -1 / 4 , possibly with logarithmic corrections, 
(iii) If initially cf> falls off as exp(— X* R x), then the associ- 
ated KPZ interface remains straight towards £ = oo. For 
this case the KPZ scaling predicts A ~ t 1 / 3 . Presum- 
ably a crossover between exponent 1/3 and 1/4 could be 
present when starting with an initial condition slightly 
faster decaying than exp(— X* R x). The identification of 
such a crossover and the modification of the global ex- [18] 
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